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Recent developments and applications of approximate actions for full lattice QCD are described. We present 
first results based on the stochastic estimation of the fermion determinant on 12'^ x 24 configurations at (5 — 5.2. 



1. BASIC PRINCIPLES 

Consider actions Si\U] and 5*2 [17] describing 
two lattice gauge theories with the same gauge 
configuration space {U} so that [i — 1, 2) 
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For example, might be the quenched Wilson 
action and 5*2 the SW-improved action for 2- 
flavour QCD. Here, F is some operator. Expecta- 
tion values in the two theories can be related via 
a cumulant expansion whose leading behaviour 
implies [|l[ 

< F >2=< F >i + < FAi2 >i +... (2) 
where A12 = 5"! - S'2, F = F- < F > etc. (3) 

In general, an action is a function of several 
parameters, e.g. the Wilson action depends on 
the bare parameters (3 and k. One may consider 
matching different actions in one of several ways 

Ml: minimise the 'distance' between the actions, 
i.e. o-2(Ai2); 

M2: match a given set of operators, i.e. require 

< Fn >l—< Fn >2] 

M3: maximise the acceptance in an exact algo- 
rithm for S2 constructed via accept/reject 
applied to configurations generated with ac- 
tion Si. 



It turns out that, to lowest order, these 3 tuning 
prescriptions coincide. They differ in a calculable 
way at next order. Details are in ||l|. 

2. APPLICATIONS TO QCD 

2.1. Tuning action parameters 

In ||l| , we demonstrated how a relatively mod- 
est number of Wilson loop operators (three or 
four) 



E 



a,W,, 



(4) 



can be used to approximate the full Wilson action 

S2^Sg-T where T = ^TrLn[MtM] . (5) 

A particularly simple application is to study 
the quenched approximation (see also |^). Here 
n — 1 and Wi ~ Sq, the usual Wilson gauge 
action. Tuning the quenched coupling /3q in Si 
to match S2 (full QCD) at a given (3 and k yields 
(at first order): 



< TWi > 

< > 



(6) 



In general, one requires a better range of loop 
operators [Q or, indeed, other types of operators 
to capture the essential physics. In all cases, the 
basic principles and techniques are the same. 
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2.2. Bare parameter dependence in QCD 

Another practical application is to take 

Si = SQCD{f3o,Ko) and S2 ^ SqcdW, k) (7) 

so as to explore the bare parameter dependence of 
the lattice theory using configurations generated 
at a finite number of reference points {Po, 1^0) in 
parameter space. For example, at fixed (3 — Po, 
one might wish to explore the k— dependence of 
measurements. According to eqn. ^ one requires 
measurements of 

A12 = Si^S2^ T{k) - T{ko) . (8) 

One can make stochastic estimates of this quan- 
tity quite efhciently (see section ^) - certainly 
much more rapidly than by performing indepen- 
dent dynamical fermion simulations at each set of 
parameters. 

2.3. Exact algorithm for full QCD 

In we constructed an algorithm which de- 
livers configurations correctly distributed with re- 
spect to a desired action 5*2 as follows: 

W2iU,U')=WiiU,U')AiU',Ai2[U']). (9) 

The acceptance probability A depends directly on 
the quality of the approximate action : 

A(Ai2)-erfc(iv/'^2(A,2)). (10) 

The action difference A12 is of course an extensive 
quantity. In a typical application, the transition 
probability Wi for the approximate action should 
correspond to a relatively fast update scheme. 
Application of this idea to full QCD may offer 
a partial solution to decorrelation problems with 
standard dynamical fermion algorithms. 

3. STOCHASTIC ESTIMATES OF THE 
FERMION DETERMINANT 

We require an unbiased estimator for TrLnH 
where H = M^M is a hermitian positive-definite 
matrix. Bai, Fahey and Golub 1^ have recently 
proposed estimators, with bounds, for quantities 
of the form 

u^f{H)v (11) 



where, for example in our case, / = Ln. Taking 
u = V as some normalised noise vector (e.g. Z2), 
we can obtain a stochastic estimate of TrLnTJ. 
The efficiency of the method, in comparison with 
the Chebychev-based methods used previously 
|^,|lj results from an elegant relationship between 
the nodes/ weights required for a A^-point Gaus- 
sian quadrature and the eigenvalues / eigenvectors 
of a Lanczos matrix of dimension N §. This re- 
lationship and resulting accuracy appears to re- 
main good, when orthogonality is lost in standard 
numerical Lanczos methods. In the present case, 
for a fixed noise vector we obtain 6 figure con- 
vergence of the quadrature with 70 Lanczos steps 
on a matrix with condition number of order 10*. 
The estimator of T = TrLni/ is of the form 

ET = —Y.I{^i), i{4,.) = Y,^Mi>^j) (12) 

i=l j=l 

where {Xj} are eigenvalues of the tridiagonal 
Lanczos matrix arising using (jji as a starting vec- 
tor and {ujj} are related to the corresponding 
eigenvectors. Details of the method will be pre- 
sented elsewhere Q. 

4. QCD RESULTS AT /3 = 5.2 

4.1. Quenched comparisons 

In Table |l] we show the results of matching 
quenched and SW-improved jsj fermion actions 
using cqns. ^, ^ The errors on comparison val- 

Table 1 

Equivalent quenched (3 results predicted from full 
QCD data at (3 = 5.2. (3q is the equivalent cou- 
pling, < P > is the measured average plaquette, 
and < P >Q the actual quenched value at [3q. 



K 


Configs. 


<P> 


I3q 


<P>Q 


.136 


100 


.4874(1) 


5.48(2) 


.488(4) 


.139 


40 


.5160(2) 


5.62(2) 


.530(4) 



ues of < P >Q include those associated with the 
uncertainty in estimating (3q . The comparison at 
the heavier quark mass, n = .136, shows com- 
plete consistency (w 0.488). That at the lighter 
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quark mass is close but not precise. It is seen from 
this and from other results that, at k = .139, one 
may be starting to see non-trivial effects of un- 
quenching. 

4.2. Static potential measurements 

To probe dynamical quark effects more closely, 
we have studied the static potential on quenched 
(/3q = 5.7) and on 2 flavour dynamical fermions 
configurations {/3 — 5.2) at several values of k. 
Measurements were made using the optimised 
techniques of Results for V{R) and the ex- 
traction of the lattice spacing via tq are pre- 
sented elsewhere |8|. Strong scale dependence on 
K is observed. We used eqns. || and || (working to 
lowest order) to 'predict' the k = .139 potential 
from that measured at k = .136 where the quarks 
are effectively much heavier. There was evidence 
that the lowest order correction captures much 
of the change in behaviour but with large errors. 
This is not surprising since making comparisons 
at fixed /3 involves different physical scales. 

4.3. Estimating 'distances' 

In Table |[ we show the distance (squared) 
cr^(Ai2), as defined in section]^, between our dy- 
namical fermion reference action at /3o = 5.2, 
kq = .136 and other relevant actions. For ex- 
Table 2 



Estimated distance squared a^(Ai2). 



Action 


P 


K 


'Distance' 


Quen. 


5.2 




9250(1330) 


Quen. 


5.479(15) 




855(255) 


Nf = 2 


5.2 


.139 


147(20) 


Nf = 2 


5.165(1) 


.139 


17(3) 



ample, the second line shows the quality of the 
quenched tuning described above (see Table |l]) : a 
reduction from 9250 to 855 by shifting /3. Also, 
we note that k = .139 is some distance (147) from 
the reference = .136. This is the origin of the 
difficulty (noted above) in predicting k = .139 re- 
sults from .136 configurations without also tun- 
ing (3. The last line of the table shows that 
operator matching should be much closer using 



(/3o, no) = (5.2, 0.136) and (/3, k) = (5.165, 0.139), 
since the estimated distance is only 17 ± 3. This 
could be checked directly with further simulations 
at /3 = 5.165, K = .139. We would expect the lat- 
tice spacing to be similar at these two points in 
action parameter space. 

4.4. Approximate loops actions 

We have begun a systematic study of possi- 
ble approximate actions constructed from Wilson 
loops. Initially, we have chosen to construct these 
from a range of link steps of size 1, 2, 3, 4 etc. The 
aim is to minimise the distance between the target 
action, e.g. full QCD at /3o = 5.2 and kq — .136 in 
the present example, and the parametrised loop 
action ||l|. It is not hard to reduce the distance 
well below that achieved with a single plaquette 
(see Table ||). However, eqn. |l^ shows, for eam- 
ple, that cr^ must be less than 10 to give an ac- 
ceptance greater than 2.5%, in the corresponding 
exact algorithm. 

Further work will include: continuing studies to 
establish more efficient loop operators and other 
classes of approximate action; studies of the pa- 
rameter space of standard dynamical fermion ac- 
tions; mapping the parameter spaces of different 
types of action. 
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